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ABSTRACT 

ON . 
Os . 

[ — . We describe a near-infrared imaging survey of Globule 2 in the Coalsack. This Bok 

globule is the highest density region of this southern hemisphere molecular cloud and 
is the most likely location for young stars in this complex. The survey is complete for 
K < 14.0, H < 14.5, and J < 15.5, several magnitudes more sensitive than previous 

Q-f observations of this globule. From the large number of background stars, we derive 

an accurate near-infrared extinction law for the cloud. Our result, Ej_h/Eh-k 
2.08 ± 0.03, is significantly steeper than results for other southern clouds. We use the 
J — H/H — K color-color diagram to identify two potential young stars with K < 14.0 

J> . in the region. We apply ff-band star counts to derive the density profile of the Coalsack 

• i-H . 

Globule 2 and use a polytropic model to describe the internal structure of this small 
cloud. For a gas temperature T ~ 15 K, this globule is moderately unstable. 

Subject headings: ISM: individual (Coalsack) — ISM: dust, extinction — ISM: globules 
— star: formation — star: pre-main sequence 



1. Introduction 

The Coalsack is a conspicuous, nearby (d ~ 180 pc; Franco 1989) dark cloud located close to 
the Galactic plane in the southern Milky Way (1 ~ 303°, b ~ 0°). Nyman et al. (1989) estimated 
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a total mass of ~ 3550 M from 12 CO data covering an area of ~ 15 deg 2 in the direction of this 
complex. Star counts toward the same area indicate an average optical extinction of Ay ~ 5 mag 
(Gregorio-Hetem et al. 1988; Cambresy 1999). However, the extinction is Ay ~ 20 mag or larger 
in several small dark globules (Tapia 1973; Hartley et al. 1986; Bourke et al. 1995a). This 
structure is similar to other well-known clouds such as Taurus, Chamaeleon I, and Lupus that are 
actively forming stars (Cernicharo & Bachiller 1984; Boulanger et al. 1998; Arce & Goodman 
1999; Cambresy 1999). Despite this similarity, previous searches of the Coalsack have failed to 
detect signs of recent star-formation (Weaver 1974; Schwartz 1977; Reipurth 1981; Nyman et al. 
1989). 

Recently Kato et al. (1999) have observed the Coalsack in the 13 CO and C 18 J=l-0 lines. 
The 13 CO data show a massive cloudlet (~ 200 M Q ) located at the western edge of the complex, 
coinciding with the region of the largest optical extinction (see Cambresy 1999). In addition, these 
authors identified five C 18 cores (see also Vilas-Boas et al. 1994), corresponding to the positions 
of well-known optically dark globules (Tapia 1973; Hartley et al. 1986; Bourke et al. 1995a). 
These cores have high column densities, typical of star-forming cores, and masses between 4-10 
M Q . However, there are no IRAS sources, and thus no evidence of recent star-formation, in these 
cores (Kato et al. 1999). 

Tapia's Globule 2 (a = 12 h 31. 5 m , S = -63° 44.5'; 2000.0) is the densest (n(H 2 ) ~ 4 x 10 3 
cm -3 ) and most massive (~ 10 M Q ) of the Kato et al. cores (see also Bourke et al. 1995b). It 
is an obvious, roughly circular patch (~ 6' radius; Bourke et al. 1995a) of extinction on shallow 
optical and near-infrared (near-IR) surveys. Jones et al. (1980) observed an area of ~ 850 arcmin 2 
to K = 9.5; Jones et al. (1984) scanned an area of 2' x 2', 70% complete to K = 13.7, with 
the central 1' x 1' region 70% complete to K = 14.7. Neither of these surveys detected candidate 
near-IR excess stars. At K = 9.5, Jones et al. (1980, 1984) could detect ~ 1 M main sequence 
stars at the distance of the Coalsack. In the central 1' x 1' region, 70% complete to K = 14.7, this 
detection limit is ~ 0.4 M Q , assuming Ak ~ 2 (see Henry &: McCarthy 1993; Delfosse et al. 2000). 

The Coalsack Globule 2 has also been observed by the Midcourse Space Experiment (MSX; 
Price et al. 2001). Fits format images and a source catalog in four mid-infrared bands, A(8.28 /im), 
C(12.13 lira), D(14.65 urn), and E(21.34 ^m), are available through the NASA/IPAC IR Science 
Archive (IRSA) 3 . Nine mid-infrared sources lie within ~ 11' from the optical center of the globule. 
One of these is located ~ 2.8' SW from the center of the most massive Kato et al. C 18 cores. 
This different is comparable to the 2' resolution of the molecular line maps. 

Here, we describe a near-IR imaging survey of ~ 15'xl5' centered on Globule 2. The obser- 
vations are uniformly complete to K = 14.0, H = 14.5, and J = 15.5. If Ak ~ 0-2 mag, these 
data are sensitive to ~ 0.2-0.5 M Q main sequence stars (see Henry &: McCarthy 1993; Delfosse et 
al. 2000). We detect ~ 6500 sources to the detection limits of K = 16.5, H = 17.0, and J = 18.0. 



3 http://irsa.ipac. caltech.edu. 
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We restrict our analysis to K = 14.0 (~ 2500 sources) where the photometric errors are relatively 
small (< 0.08 mag). 

In §2 we describe the observations and data reduction. In §3 we derive a reliable extinction 
curve for background stars in the Globule 2 region and use this result in §4 to search for potential 
young stellar objects in our survey region based on their locations in the J — H/H — K diagram. 
We identify only two objects with K < 14.0 and near-IR excesses characteristic of young stellar 
objects. In §5 we use If-band stars counts to derive the density profile of this small cloud and 
compare it with the predictions of polytropic models. Assuming T ~ 15 K, Globule 2 is moderately 
unstable. We conclude with a brief summary in §6. 

2. Observations and Data Analysis 

We obtained JHK imaging data for Globule 2 in the Coalsack and four relatively unreddened 
control fields on 13-16 February 1995, and 8-11 March 1996, with CIRIM (the Cerro Tololo Infrared 
Imager) at the CTIO 1.5m telescope. The CIRIM uses a 256 x 256 HgCdTe NICMOS 3 array, 
which provides a field of ~ 4.9' x 4.9' with a plate scale of 1.16" per pixel. We covered an area of 
~ 15' x 15' on the cloud on a regular grid, with 1' overlap between adjacent frames. We acquired 
two 6x5 sec exposures for each field, shifted by 20". We also observed four 5' x 5' control fields, 
selected from a visual inspection of the ESO Red Sky Survey prints, close to Globule 2 region and 
relatively free from significant optical extinction. Table 1 lists positions of these off-cloud regions. 

We process the data using standard techniques with the IRAF 4 software package. Gomez Sz 
Kenyon (2001) describes the data reduction procedure in detail. Briefly, we apply dark, flatfield, 
and sky corrections to all images, align dithered pairs using the IRAF subroutines GEOMAP and 
GEOTRAN, and then add the dithered pairs. We select the K images as the reference frames and 
transform the J and H images to the same pixel scale as the K images. We use DAOFIND to 
locate stars 4a above the local background and add to the DAOFIND list all stellar objects missed 
by this routine and found by visual inspection of each image. We then derive photometry for each 
image using the APPHOT PHOT task, using a circular aperture with 5" radius. We detect ~ 6500 
sources at JHK in our survey region. 

To calibrate our photometry, we observed on each night a set of 10-15 standards from Elias et 
al. (1982) and from the UKIRT faint JHK standard stars list (Hawarden et al. 2001). Gomez & 
Kenyon (2001) estimate photometric uncertainties of ± 0.03 mag in this calibration. These authors 
also list typical photometric uncertainties in each magnitude bin. These uncertainties reach 0.3 mag 
for J ~ 17, H ~ 16, and K ~ 15. The 5a limiting magnitudes are K = 14.5, H = 15.5, and J = 
16.0. The magnitudes are in the CIT (Elias et al. 1982) system (see Gomez & Kenyon 2001). 
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To check the completeness of the survey, we construct a list of detected sources. Figure 1 shows 
the magnitude histograms of background stars as functions of JHK. The number of stars per bin 
increases monotonically up to the completeness limit and then turns over. The corresponding limits 
in each filter are J = 15.5, H = 14.5, and K = 14.0. 

To derive coordinates for our survey stars, we use WCSTools 5 (Mink 1997), a suite of pro- 
grams to calculate a direct transformation between the coordinates of the image (x,y) and the sky 
coordinates (a, 6). We measure the transformation coefficients adopting matches between program 
stars and stars in the U.S. Naval Observatory SAl.O Catalogue (Monet et al. 1996). We estimate 
an average uncertainty in our positions of 1". This estimate relies on a direct comparison of our 
coordinates with coordinates from the Digitized Sky Survey 6 (DSS). 

3. The Coalsack Globule 2 Near-Infrared Reddening Law 

We use a general photometric technique to derive the near-IR reddening law in the Coalsack 
(Kenyon et al. 1998). The technique assumes that the stellar population behind the cloud is 
identical to the stellar population in off-fields several degrees away. Because the Coalsack contains 
so few pre-main sequence stars, this assumption is reasonable. Kenyon et al. (1998) derive the 
J — H and H — K color excesses for each on-cloud source relative to every off-cloud source and then 
compute the average and median color excesses for each on-cloud source. The probable error of 
the average color excess is the sum in quadrature of the errors of the on-cloud and off-cloud colors. 
For the median color excess, the probable error is the inter-quartile range. 

Figure 2 shows average color excesses of Coalsack stars with K < 14. We consider stars with 
K < 14.0, where the photometric errors are < 0.1 mag. We eliminate two potential young stars 
with near-IR excesses from this analysis (see §4 below). The color excesses are correlated, with a 
Spearman rank-order correlation coefficient of r s = 0.88. The probability for no correlation between 
the two color excesses is formally zero according to the Spearman rank-order test. Straight line fits 
to color excess measurements yield Ej_h/Eh-k = 2.07 ± 0.05. Median color excesses for each 
source yield the same slope and error. 

To test the robustness of the near-IR extinction law, Kenyon et al. (1998) define the reddening 
probability function. The reddening probability p(E x , E y ) is the chance of measuring a pair of color 
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Fig. 1. — Histogram distribution of magnitudes for background stars in our survey region. The 
completeness limit in each band is indicated by the dashed vertical lines. 
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Fig. 2. — Average color-excess diagram for Globule 2 stars with K < 14. The sample includes 2271 
stars used in the analysis. 
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excesses E x and E y , where x and y are color indices and N = f p(E x , E y ) dE x dE y is the number of 
reddening measurements. The number, N, is also the number of on-held stars in the sample. The 
density function, p(E x ,E y ), depends on the distributions of colors in the on-held and the off-held 
and results from a convolution of the densities of off-held and on-held stars in the color-magnitude 
diagram (Kenyon et al. 1998). These color distributions use a kernel density estimator with 
smoothing parameter h = 0.2 to derive the density. 

Figure 3 shows the probability functions for the Coalsack. The color-color diagram in the left 
panel shows the middle of the main sequence and the lower main sequence as peaks in the density. 
The middle panel indicates a clear peak in the sample of Coalsack stars at H — K f» 0.5 and 
J — H Ki 1.3-1.5. The reddening probability function in the right panel consists of a collection of 
nested aligned ellipses. The major axes of these ellipses have a slope consistent with the slope of 
the reddening law in Figure 2, Ej_h/Eh~k = 2.08 ± 0.03. Kenyon et al. (1998) describes the 
derivation of the slope of the major axis of each contour in Figure 3. The ellipses are symmetric 
about the reddening line, suggesting that the Coalsack contains few stars with near-IR excesses 
(see below). 

The Coalsack near-IR extinction law is different from other near-IR extinction laws. In previous 
studies, we derive Ej_h/Eh-k = 1-80 ± 0.03 for the Chamaeleon I dark cloud (Gomez & Kenyon 
2001) and Ej„ H /E H _ K = 1.57 ± 0.03 for the p Ophiuchi dark cloud (Kenyon et al. 1998). The 
Cha I extinction law agrees with the standard extinction law of Bessell &: Brett (1988), Ej_h / Eh~k 
= 1.75, transformed to the CIT system. The p Oph law is more similar to results from the He et 
al. (1995) survey of luminous southern stars, Ej^h/Eh~k = 1-47 ± 0.06. The extinction laws 
for the three dark clouds differ by ~ 10u, suggesting dramatic changes in grain composition or size 
between these clouds. 

The ratio of total to selective extinction, Ry , may also to vary among dark clouds (see Kenyon 
et al. 1998, and the references therein). Kim & Martin (1996) attribute this variation to different 
grain properties from cloud to cloud. Dense regions have larger Ry and fewer small grains than 
diffuse clouds. On average, grain sizes increase with cloud density. In addition, recent calculations 
show a clear dependence of the extinction law on the grain mixture and composition (see, for 
example, Vaidya et al. 2001). 

Molecular line observations confirm the difference between Cha I and the Coalsack. Harjunpaa 
& Mattila (1996) report a signihcant range in the ratio of the CO column density N(CO) to the 
J — K color excess Ej_k in several southern molecular clouds. The observed N(CO)/Ej-k is 
a factor of two larger in Cha I than in the Coalsack. The R CrA star-forming region also has a 
large N(CO)/Ej_ K . Harjunpaa & Mattila (1996) suggest that N(CO)/N(H 2 ) or N(H 2 )/Ej_ K 
or both vary signihcantly between these clouds. Other clouds display similar variations. Kramer 
et al. (1999) note changes in N{CO)/Ay within the dark cloud IC 5146, which they interpret as 
differences in the amount of CO depletion onto dust grains. From observations of C 17 0, C 34 S, and 
N2H+, Bergin et al. (2001) show that the pattern of gas depletion in IC 5146 is consistent with 
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Fig. 3. — Density distributions for Globule 2 stars. The left and middle panels show the color density 
distributions for stars off the main Globule 2 field (off-field, left panel) and for stars projected on 
the globule (Globule 2, middle panel). The right panel plots the reddening probability distribution 
for the Globule 2 sources, assuming that these stars have the same color distribution as the off-field 
stars. 
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models of cloud chemistry. 

Our observations of p Oph, Cha I, and the Coalsack indicate important differences in grain 
chemistry among dark clouds. Changes in the slope of the near-IR extinction law appear correlated 
with the degree of star formation activity or the density of the cloud or both, p Oph is the most 
active star formation region of this group and it has the regions of largest extinction. The Coalsack 
is the most quiescent cloud and has the lowest extinction. Testing whether this correlation is real 
or a fluke of small number statistics requires high quality extinction measurements in more dark 
clouds. 



4. Search for Candidate Pre-Main Sequence Stars 

Figure 4 shows the H — K / J — H diagram for our complete survey region. The solid line 
corresponds to the main sequence locus (Bessell &, Brett 1988); the dashed lines, parallel to the 
reddening vector (Ej_h/Eh-k = 2.08 ± 0.03) derived in §3, define the reddening band extending 
from the main sequence. We detect ~ 2500 stars in Globule 2 with K < 14.0. 

We distinguish two groups of sources in Figure 4: i) stars that lie along or follow the reddening 
band direction, and ii) two sources with near-IR excess located to the right of the reddening 
band. The first group contains almost all of our detections. These objects are reddened main 
sequence or giant stars (see §3) and lie behind the cloud 7 . The color excesses of the second sample 
cannot be attributed to reddening by normal dust in the cloud. Circumstellar material (disks + 
infalling envelopes) surrounding the young central stars can produce these color excesses at near-IR 
wavelengths (Lada & Adams 1992; Kenyon & Hartmann 1995). Two stars in our sample have near- 
IR color excesses similar to previously known young stars Table 2 lists coordinates and magnitudes 
for these stars. They are located within ~ 5-6' of the position of the C 18 core detected by Kato 
et al. (1999). In a similar manner, the nearest mid- infrared MSX sources to these candidates lie 
at distances of ~ 5' and 3.5', respectively. 

We find three objects in common between our near-infrared survey and the nine sources from 
the MSX catalog that lie within radius of ~ 11' from the optical position of the Globule 2. The 
other 6 mid-infrared sources are located close but outside of our 15' x 15' observed area. 

Although the Coalsack lies in the general direction of the Sagittarius-Carina arm of the Milky 
Way, Nyman et al. (1989) note that there is no background cloud at the position of Globule 2. 
Confusion with background sources is not significant; the selected candidates probably belong to 
the Coalsack region lying in the vicinity of this globule. The near-IR excesses of these sources 
suggests the pre-main sequence nature of these stars. Optical spectroscopy is required to test this 
suggestion. 



Any weak emission T Tauri stars that belong to the cloud and that, as a class of young stellar objects, have no 
significant near-IR excess are included in this group. 
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Fig. 4. — Color-color diagram for near-IR sources with K < 14.0 in our survey region. The solid 
line indicates the main sequence locus (Bessell & Brett 1988) and the two parallel lines define the 
reddening band, extending from the main sequence. To define this band we use the reddening 
vector, Ej-h/Eh-k = 2.08 ± 0.03, derived in section 3. Typical photometric errors for K ~ 
14.0 are displayed in the upper left corner. The big stars show the positions of two near-infrared 
candidate young stellar objects. 
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5. The Coalsack Globule 2 Density Prone 

5.1. The Extinction Map 

To map the distribution of extinction in Globule 2, we divide the observed region (a 15' x 
15' area) into a fixed grid of 1' x 1' overlapping squares, separated by the 30" Nyquist spatial 
sampling interval. We follow the same procedure for our off-cloud regions (see Table 1). ii-band 
star counts in these two areas then yield an extinction map. This standard method compares the 
stellar density in each position of the cloud with the density in a nearby control area assumed free 
of significant obscuration. The number of stars per unit area brighter than in the control field, 
N ff, is 



Log N off = a + b m A . (1) 

On the cloud, the slope of equation (1) remains the same but the number of stars per unit 
area, N on , decreases as the extinction Am\ increases. The extinction is then 

A 1 T N off ,„v 

= b L ° g Non ' (2) 

where b is the slope of the cumulative luminosity function. We use the ff-band luminosity function 
for off cloud regions to measure b = 0.32 ± 0.01. We convert the -ff-band extinction, Ah, to optical 
extinction, Ay, using the Rieke & Lebofsky (1985) reddening law (Ay = 5.7 Ajj) and convolve 
the individual stellar measurements with a two dimensional Gaussian filter with a FWHM of 60", 
the size of our grid cells. Large cells yield poor resolution; smaller cells have too few stars for 
a meaningful measurement. Figure 5 shows our extinction map. The smoothed stellar positions 
display roundish shapes. The outermost contour corresponds to Ay = 5. Subsequent contours 
increase in steps of 3 mag, reaching Ay = 23 mag in the innermost area. As noticed by Jones et 
al. (1984), Globule 2 is embedded in the Coalsack; the 5 mag level is probably the average level 
of dust absorption in the surrounding cloud region. 

We determine mean radii for each optical obscuration level, normalize this 'average' density 
with respect to the 23 mag peak level, and plot the results in Figure 6 as dot symbols. The projected 
radial dust column density distribution is shallow, suggesting a configuration near equilibrium. 
Density profiles of unstable clouds show a clear enhancement toward the center, indicating the 
collapse of material (see Harvey et al. 2001). 



5.2. The Polytropic and the Bonnor-Ebert Models 



To model the internal structure of Globule 2, we consider the Lane-Emden equation for a 
self-gravitating isothermal gas sphere in hydrostatic equilibrium: 
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Fig. 5. — Extinction map for the Globule 2. The lowest contour corresponds to Ay = 5 mag. 
Subsequent contours increase in step of 3 mag. The coordinates of the center of the region are: 
a(2000.0) = 12 h 31 m 24 sec , 5(2000.0) = - 63° 45' 11". The big cross indicates the position of the 



-13- 



0.8 



T3 

S 0.6 



to 

S 0.4 



0.2 



Polytropic Model (n = 2.3) 



0.2 



0.4 0.6 
Fractional Radius 



0.8 



0.8 



t3 0.6 

CD 



0.4 
0.2 



Bonnor-Ebert Model (p = 2.5x10 g cm ) 



15 K — L 



7.0 ± 0.3 



0.2 



0.4 0.6 
Fractional Radius 



0.8 



Fig. 6. — Normalized density profile for the Globule 2. The optical extinction Ay is referred to the 
central peak value (23 mag). The radial distance is the fractional radius with respect to R = 140" 
(~ 0.1 pc), corresponding to Ay = 5 mag. The dot symbols represent our observational data and 
the continuous lines the model fits. The upper panel corresponds to a polytrope with n = 2.3. 
The lower panel indicates the Bonnor-Ebert model for a central density p c = 2.5 x 10 -19 g/cm 3 , 
assuming T = 15 K for this globule. 
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where £ = r/a and 9 = T/T c are the Emden variables, with a an unspecified factor having 
dimensions of length and T c the central temperature (see, for example, Novotny 1973). The 
parameter 'n' is the polytropic index relating the gas pressure to the density, P gas = K n p n . We 
consider two approaches to solve equation (3): a) a numerical integration for specific values of 'n' 
and b) the introduction of a variant of this pressure-confined, self-gravitating isothermal sphere 
equation. 

Kenyon & Starrfield (1979) adopted the first approach and numerically solved equation (3) for 
different values of the polytropic index n, satisfying the following boundary conditions: 9=1 and 
d#/d£ = at ^ = 0. The upper panel of Figure 6 shows the normalized density profile superposed on 
a polytropic model with index n = 2.3. We adopt an 'average' external radius of 140", 0.1 pc at the 
distance of the this cloud. This radius corresponds to Ay = 5 mag, which we assume represents the 
average optical extinction in the surrounding Coalsack region. For n = 2.3., Kenyon & Starrfield 
(1979) 's calculations predict a central density p c = 3 x 10~ 19 g/cm 3 and a temperature T = 20 
K to support the cloud against collapse. Although we have no direct temperature measurement of 
Globule 2, the average temperature for other small globules is T = 10-15 K (Bourke et al. 1995b). 
The factor of 1.2-2.0 difference between the 'observed' and hydrostatic equilibrium temperature 
suggests that Globule 2 is marginally unstable and on the verge of collapse. 

In our second approach, we consider a variant of equation (3), known as the Bonnor-Ebert 
model (Bonnor 1956; Ebert 1955) or the modified Lane-Emden equation: 



il( f * = e -* (4) 

ed^ de j ' w 

where p = p c e~^\ Here, £ = (r /a)^4irGp c is a dimensionless radial parameter and a is the 
isothermal sound speed (a = ^kT/pm-n; p is the mean molecular weight and me is the mass of 
the H atom). Introducing standard boundary conditions (0(0) = and d^(0)/d£ = 0), we solve 
equation (4) numerically. Assuming an external radius R (140" or 0.1 pc), solutions of equation 
(4) can be parametrized by p c . The lower panel of Figure 6 shows our best fit for p c = 2.5 x 10~ 19 
g/cm 3 superposed on the Globule 2 density profile. This solution assumes T= 15 K as a typical 
temperature for Bok globules (see Bourke et al. 1995b) and p = 2. 

Bonnor (1956) demonstrated that the value of £ at the external radius (R), £ max = R/a \J AirGp c 
provides a stability measure. Systems with £ max > 6.5 are unstable to gravitational collapse. For 
Globule 2, £ max = 7.0 ± 0.3 exceeds the critical value and is moderately unstable. 



The mass of the globule and the pressure at its outer boundary are functions of £: 



(5) 



P(0=4e"^, (6) 

where f3 = a/ ^/4:TrGp c . Evaluating these expressions at r = R (or £ = £ max ), we derive M = 4.5 M Q 
and P cxt = 9 x 1CT 12 dyn/cm 2 . Kato et al. (1999) estimate 10 M ; our external pressure is 
roughly about 6.5 times the value for the interstellar medium (see Elmegreen 1991). 

Both the polytropic and the Bonnor-Ebert models predict a moderately unstable configuration 
for Globule 2. The 'typical' T ~ 10-15 K is a factor of 2 lower than the minimum value required 
for hydrostatic equilibrium. Magnetic fields may help to prevent the collapse. Jones et al. (1984) 
found evidence of a mildly compressed magnetic field from .ff-band polarization measurements of 
background stars in the Globule 2 region. The very modest intensities of magnetic fields associated 
with dark clouds (B ~ 16 pG; Crutcher et al. 1993) suggest, however, that magnetic fields can 
not prevent the collapse indefinitely and alter the global process significantly. 

Alves et al. (2001) used the Lane-Emden equation to fit the density profile of Barnard 68 
(Bok 1977), an isolated nearby (d = 125 pc) southern hemisphere Bok globule. They derived 
£max = 6.9 ± 0.2, slightly larger than the critical value. B68 has a radius of 0.06 pc, M = 2.1M Q , 
T = 16 K, and an external pressure of 2.5 x 10~ n dyn/cm 2 . This starless globule shows no 
evidence of star-formation and thus may be about to collapse to form stars. 

Harvey et al. (2001) analyzed the dense globule in B335 (d ~ 250 pc; Tomita et al. 1979). 
The shape of the B335 density profile is steeper than either Globule 2 or B68. Harvey et al. (2001) 
show that collapse models and an unstable Bonnor-Ebert sphere (with £ max = 12.5 ± 2.6) match 
the observed density profile. Molecular line profiles provide additional evidence that this globule 
is collapsing (Zhou et al. 1993; Choi et al. 1995). B335 is associated with a deeply embedded 
young stellar object (Keene et al. 1983) that drives a bipolar outflow (Cabrit et al. 1988; Hirano 
et al. 1988, 1992). Thus this core is in a more advanced evolutionary stage than B68. 

The Coalsack Globule 2 may represent an intermediate phase between B68 (a pre-collapse 
globule) and B335 (an already collapsing core) according to the polytropic models and Bonnor's 
critical parameter, £ max . 

6. Summary 

Our near-IR survey of Globule 2, the highest density and more massive core in the Coalsack 
complex, leads to three primary conclusions. 

The slope of the near-IR extinction law in the cloud, Ej-h/Eh-k = 2.08 ± 0.03, is much 
steeper than in p Ophiuchi or Chamaeleon I. If real, the trend of increasing slope with decreasing 
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star-formation activity in these three clouds suggests changes in grain chemistry with extinction or 
star formation activity. 

Based on near-IR excess emission, we detect two pre-main sequence candidates with K < 14 in 
the vicinity Globule 2. If confirmed as pre-main sequence stars with optical or near-IR spectroscopy, 
these are the first pre-main sequence stars discovered in the region. The globule contains no known 
IRAS sources (Nyman et al. 1989; Bourke et al. 1995a) or other pre-main sequence candidates. 
Our candidates are too faint for detection with IRAS. The low success rate in finding pre-main 
sequence stars in the Globule 2 region illustrates the low activity of the Coalsack as stellar nursery, 
in agreement with previous investigations (Weaver 1974; Schwartz 1977; Reipurth 1981; Nyman 
ct al. 1989). 

We use -ff-band star counts to derive the density profile of this Bok globule. For a 'typical' gas 
temperature T ~ 15 K, model fits suggest this small cloud is moderately unstable, with a Bonnor 
critical parameter £ max = 7.0 ± 0.3. The mass derived from these models, M = 4.5M Q , agrees 
with estimates derived from the CO column density. 

We are grateful to the CTIO staff, specially to M. Fernandez, M. Hernandez, and P. Ugarte for 
assistance during the observing runs. We also thank R. Elston and R. Probst for their help with 
CIRIM and D. Mink for assistance with the WCSTools software. An anonymous referee carefully 
read the original manuscript and provided useful suggestions. 
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Table 1. Off-Cloud Regions. 



Region 


a(2000.0) 


(5(2000.0) 


I 


12 32 53 


-64 21 34 


II 


12 32 55 


-64 27 43 


III 


12 33 46 


-64 26 34 


IV 


12 33 49 


-64 22 44 



Note. — Units of right ascension 
are hours, minutes, and seconds, 
and units of declination are degrees, 
arcminutes, and arcseconds. 



Table 2. Candidate Pre-Main Sequence Stars with K < 14.0. 



Star 


a(2000.0) 


5(2000.0) 


K 


H -K 


J -H 


1 


12 31 54.4 


-63 50 00 


13.99 


0.66 


0.90 


2 


12 32 05.0 


-63 49 30 


12.56 


0.45 


0.72 



Note. — Units of right ascension are hours, minutes, and 
seconds, and units of declination are degrees, arcminutes, 
and arcseconds. 



